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. . . Abstract 

P^ ' Relaxation and decoherence of a qubit coupled to environment and driven by a resonant ac field are investigated by analytically 

^^ , solving Bloch equation of the qubit. It is found that the decoherence of a driven qubit can be decomposed into intrinsic and 

C^ ■ field-dependent ones. The intrinsic decoherence time equals to the decoherence time of the qubit in free decay while the field- 
dependent decoherence time is identical with the relaxation time of the qubit in driven oscillation. Analytical expressions of 

^..^, the relaxation and decoherence times are derived and applied to study a microwave-driven SQUID fiux qubit. The results are 

C^ ■ in excellent agreement with those obtained by numerically solving the master equation. The relations between the relaxation 

^ ■ and decoherence times of a qubit in free decay and driven oscillation can be used to extract the decoherence and thus dephasing 

—< J times of the qubit by measuring its population evolution in free decay and resonantly driven oscillation. 
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I. INTRODUCTION 

In recent years, significant progress has been made on physical iniplcmcntation of quantum computation based 
on superconducting qubits. Quantum coherence has been successfully demonstrated in a variety of superconducting 
single-qubit systems [ll, 0, S 0, B H S S] and coupled two-qubit systems, [§, [lO, \lM, Il4| indicating the potential 
of superconducting qubits in quantum computing. However, due to unavoidable coupling with environment, the 
superconducting (charge, flux, and phase) qubits always suffer from decoherence such as relaxation and dephasing, 
resulting in relatively short coherence times. For this reason, environment-induced decoherence has been and is still 
a main obstacle to the practical application of superconducting qubits in quantum computation. [2, H^ [ij, Ha \lM 

The environment-induced decoherence of superconducting qubits has been extensively studied both theoretically 
ISJI^'JI^'M E i20, 21, 21, 23, 24, 25, 2Q, 21, 28, 29, 30, 31, 32, 3| and experimentally @, i, JH, Ml, 25, 34, 35, M 
4011 i n the absence of ac driving fields (free decay). Quite a few proposals, such as dynamical decoupling, 
H, [45l. l46l - 1471 l48j decoherence free subspaces, [43, S, \EM, 113 spin echoes, [3, 0, [SJI and coherence- 
preserving qubits, [53J have been proposed to reduce such kind of decoherence. However, in superconducting-qubit 
based quantum computation ac fields (e.g., microwave fields) are usually used to manipulate the qubit's state. 
[H, 0, S 0, S B 0, @, M, [l3, \1M, E, [3 Recent experiment [ij] shows that the decoherence time of a superconducting 
qubit is significantly increased in the presence of a resonantly ac driving field. Thus a comprehensive understanding 
to decoherence of a realistic superconducting qubit needs to include influence of driving fields. [iSl 

In this paper, we study the effect of driving fields on the relaxation and decoherence of a driven qubit. We 
focus our study on weak resonant driving fields characterized by fJ <C w^^ (where, fi and w^ are the Rabi frequency 
and the frequency of ac driving field, respectively) and low temperatures. In this case, population leakage to non- 
computational states due to strong field effect |54l . |55| and thermal activation 56] is negligible and thus a multilevel 
superconducting qubit coupled to environment and driven by a resonant ac field can be well approximated by a 
resonantly driven dissipative two-level system (TLS). We first explore relaxation and decoherence of a qubit in free 
decay and demonstrate that our results agree very well with those in the literature. We then derive analytical 
expressions of relaxation and decoherence times for a driven qubit through analytical solutions of Bloch equation 
of the driven dissipative TLS. The relations between the relaxation and decoherence times can be used to extract 
decoherence and dephasing times of the qubit by measuring its population evolution in free decay and driven oscillation. 
Finally, we use the analytical expressions to study relaxation and decoherence of a superconducting quantum interface 
device (SQUID) flux qubit driven by a microwave and show that the analytical results are in excellent agreement with 
the results obtained by numerically solving the master equation. 

II. BLOCH EQUATION OF A GENERAL DRIVEN DISSIPATIVE TWO-LEVEL SYSTEM 

A dissipative system can be described by a reduced density operator. In the Hilbert space spanned by the eigenstate 
|m) of eigenenergy £'„j (m = 1, 2, • • • ), the reduced density operator p is represented by a reduced density matrix with 
matrix elements p„„ = {m |p| n) . The diagonal matrix element Pmm and off-diagonal matrix element p„„ {m ^ n) 
are the population and coherence of the system, respectively. 

In general, the reduced density operator is governed by generalized master equation of non-Mar kovian process. 
|23l |33 . ISTI . ISSJ I However, in the case of weak damping and low temperature considered here, the generalized non- 
Markovian master equation is equivalent to the Markovian master equation |32l | and in the case of resonant driving, 
the results obtained from the Markovian master equation are the same as those from the non-Markovian master 
equation. [57.] Furthermore, in the case of weak damping. Lamb shifts are usually very small compared to the ac 
driving field and thus are often neglected [53|- After dropping the terms related to the Lamb shifts, the Markovian 
master equation for the reduced density matrix of a driven multilevel qubit is cast into the generalized Bloch- Redfield 
equation. [5a, IfiOJ, |6l| For a driven dissipative TLS with conserved population, the reduced density matrix elements 
satisfy P21 = PI2 and Pn + P22 ~ 1 ^^'^ ^^^ governed by [5y, [6i| 

dt 

+ ^11,12(^12 +P21) +-^11,22^22, (1) 



*^12 (P2I ~ P12) ^ -^22,11^11 
-Rll,12 {Pl2 + P21) + -^11,22, 

= «^2lPl2 - * [^12 (P22 - Pll) 



^Pl-2 _ „■ . „ a\tjF 



dt 

+ (-f^ll ^ ^22) P12] + ^1241/3ll 

+ ^12,12Pl2 + Rl2,2lP2l + -Rl2,22P22i (2) 

where, tOmn = {Em — En) /h is the transition frequency, Rmn,m'n> is the matrix element of damping rate superoperator 



describing the effect of environment on the TLS, [Sa, [61| and ff,^„ = {m \Hp \ n) /h is the matrix element of interaction 
Hamiltonian Hp between the TLS and resonant ac driving field. For simplicity and without loss of generality, 
henceforward, unless otherwise specified, we assume 7?/^ = ecos(w^<) and H(i — H22 — 5 cos {(jj^t), where e and 5 
are the constants associated with transition matrix elements and field strength. 

The damping of a dissipative system is caused by the interaction between the system and environment such as a 
thermal bath. In the case of weak damping, the interaction Hamiltonian Hi between the system and the thermal 
bath can be approximated by a linear function of system variable and expressed as [la. I3ll.l60l| Hj = —AxUb, where, 
X is the generalized coordinate of system, Ub is a function of bath variable(s), and A is a coupling constant. The 
effect of the thermal bath on the system can equivalently be characterized by a spectral density J (lo) . [62, mM I'^ this 
case, the matrix elements of the damping rate superoperator can be calculated by [50, mM 
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(3) 



where, Xmn is the transition matrix element. 

An equivalent approach to describe the TLS is to use Bloch vector {u,v,w) defined by 

u = TripcTx) = 2Repi2i 
V = Tr (per y ) = -2Impi2, 
w = Tr{pa^) =pii -p22- 



(4) 
(5) 
(6) 



The components of the Bloch vector u and v represent the real and imaginary parts of the coherence Pi2, while the 
component w is the population difference. Substituting the Bloch vector into the master equation ([1]) and ([2]) one 
obtains the following Bloch equation 



du 
'dt 

dv 
'dt 



dw 



= (^1242 + -Rl2,2l) U + [UJ21 - (5 cos (uj^t)] V 

+ (-Rl2,ll ~ Rl2.22)'W + (-Rl2ai + -^12, 22) , 
= - [W21 - 5 cos {Ulfj_t)] U + (i?12,12 - -Rl2,2l) V 

—2e COS {lj ^t) w , 



— = 2Rii^i2U + 2eCOs{uJ^t)v-{R22,ll+Rll,22)w 
+ (^11,22 ~ -^22,11) • 



(7) 
(8) 
(9) 



III. RELAXATION AND DECOHERENCE IN THE ABSENCE OF AC DRIVING FIELDS 



Suppose the initial value of the Bloch vector is {uo,vo,wo). In the absence of ac driving fields (free decay) and 
under the rotating- wave approximation (RWA) , the solution of the Bloch equation is given by (see Appendix A for 

details) 



u = [Acos{-cut - 9 q) + B sin {-cut)] e~'^*, 
V = [-Asiii^mt) + Bcos{mt + 6o)]e^''\ 



(10) 

(11) 
(12) 



where, k = — i?i2,i2 is introduced for convenience, 7 = i?22,ii + ^11,22 is the relaxation rate in free decay, Woo, '^, and 
^0 are the parameters related to the damping rate matrix elements, and A, B, and D are the constants determined 
by the initial condition of the TLS and damping rate matrix elements. These parameters and constants are given 
by Eqs. (|A4[) to (|A9[) in Appendix A. It is shown that the population and coherence undergo simple exponential 
decays independently at rates 7 and k, respectively, and the coherence components u and v contain fast oscillating 
factors of frequency tu. When i — > 00 the coherence components tend to zero, while the population difference tends 
to a constant Wac,. In the case of zero damping (i?i2.2i = 0) w = UJ21 and 6q — Q while in the case of weak damping 



(^12,21 "^ ijJ2i) ^ < 1^21 and 60 > 0. Hence Woo is the population difference at t — > 00, tu is the transition frequency 
modified by the damping, and Oq is the phase shift due to the damping. 

It is shown from Eq. (fT2|) that the population (difference) decays exponentially with the rate of 7. Hence, the 
relaxation time is 

Ti^j-\ (13) 

From Eqs. (fTO|l and (fTTjl the coherence, Pi2 = (" ^ iv)/2 ex e^"*, decays exponentially with the rate of k. The 
decoherence time is therefore 

T2 = K-\ (14) 

These results agree very well with those obtained by others. JSo] Eqs. ([T3|) and (fT4|) indicate that for free decay Ti is 
independent of T2. In general T2 < 2Ti due to dephasing. [33| The pure dephasing time Tip is related to Ti and r2 
by© 

111 

— = • (15) 

T^ T2 2Ti ^ ' 

For the general case of a qubit coupled to a thermal bath, using Eq. ([3]), we obtain 

^r' = ^|j ("=.) + ^("12)1, (16) 

and 

These results demonstrate that the relaxation and dephasing rates of the a qubit are determined by the s pec tral 
densities at the transition frequency uj — \uj2i\ and low frequency u; — 0, respectively. [IJ, [I^, [Sll IH, |33, Ha] In 
addition, the relaxation rate is proportional to the modulus square of the transition matrix element |a;i2| , while the 
dephasing rate is proportional to the square of the difference of average coordinates of the two states {xn — a;22) ■ 
For a qubit having {xn — X22) = 0, the dephasing is completely suppressed. 

IV. RELAXATION AND DECOHERENCE IN THE PRESENCE OF AN AC DRIVING FIELD 

For a driven dissipative TLS, e ^ and also (5 7^ in general. Under the RWA, the solution of the Bloch equation 
in the underdamped regime with e > rj is given by (see Appendix B for details) 

u = Acqs {rut - 9o) e^*** + Bo sin (rut) 

+Bism{nt + 0i)sm{mt)e^^\ (19) 

V = — A sin {zut)e^'^'^ + Eq cos {iut + 9 0) 

+Bi sin {nt + 6*1 ) cos {vut + 6*0 ) e"^* , (20) 

w := Da + Disin{nt + 92)e-'^\ (21) 



where, ry = (7 — k) /2 is the damping strength, fi = i/e^ — 77^ is the Rabi frequency, F = (7 + k) /2 is the mean value 
of relaxation and decoherence rates in free decay, A is the field-independent constant identical with that in Eqs. (fTU)) 
and (fTTj) and is given by Eq. (jA7p in Appendix A, Bq, Bi, Dq, Di, 9i, and 02 are the field-dependent constants given 
by Eqs. (jB4[) to (jBlip in Appendix B. For quantum information processing, manipulation of qubits must be in the 
underdamped regime. Hence, hereafter our discussion is focused on the resonantly driven qubit in this regime. 

It is shown clearly that the solution given by Eqs. p^ to (PT|) contains a slow Rabi oscillation. It is also shown that 
in general the population and coherence undergo more complicated damped oscillations which are totally different 
from those in free decay. The first terms of the coherence components on the right-hand sides of Eqs. P^ and PD)) 
are independent of the driving field. They only depend on the initial state and matrix elements of the damping rate 
superoperator. Hence they represent intrinsic decoherence induced by the coupling between the qubit and thermal 



bath. Note that the intrinsic decoherence rate k of a qubit in driven oscillation equals to the decoherence rate of free 
decay (see Eqs. (fTO|) and (fTTjl ). In contrast, the remaining terms of the coherence components on the right-hand sides 
of Eqs. ([T9| and ([20| strongly depend on the driving field. In particular, the third terms decay exponentially and 
thus represent field-dependent decoherence. The field-dependent decoherence rate T and the relaxation rate of the 
qubit in driven oscillation are equal (see Eq. ((2T|) ). 

Note that both the relaxation and decoherence should have been eliminated in principle ii A = Bi = Di = Q in 
Eqs. ()19|) to (I2ip . However, such conditions can not be fulfilled in reality for a resonantly driven qubit starting from a 
pure state. Thus the relaxation and decoherence can not be suppressed simultaneously with a single resonant driving 
field. [42, |4j, 0, Hj Sa, |43| Nevertheless, as will be demonstrated below, for a special class of initial states, a properly 
chosen resonant driving field could slow down the decoherence. [6J] 

Since in general the population and coherence of a driven qubit decay with more than one rates, it is necessary to 
introduce multiple relaxation and decoherence times to describe the system's behavior properly. Each relaxation time 
or decoherence time represents one of the characteristic times of the qubit for each specific case. For instance, from 
Eq. ([2T|) , the relaxation time Ti is defined by 

Ti=r-i, (22) 

where the tilde "~" is used to indicate that the qubit is resonantly driven. Similarly, from Eqs. (fT9| and ((20|) . the 
intrinsic decoherence time T2,i and field-dependent decoherence time T2.2 are defined by 

%,! = 1^-' = ^2, (23) 

and 

T2^2 = r-i = ri, (24) 

respectively. Namely, the intrinsic decoherence time T2,i and the field-dependent decoherence time T2,2 of a qubit 
in driven oscillation equal to the decoherence time T2 of the qubit in free decay and the relaxation time Ti of the 
qubit in driven oscillation, respectively. If the initial state of the driven qubit corresponds to uq = (e.g., the ground 
state), then the intrinsic decoherence terms vanish and the decoherence of the qubit is characterized completely by 
the field-dependent decoherence time T2.2- Note that the Eqs. ([14]) and (|22p are used to obtain the second equality 
of Eqs. ([23|) and ([24|). 

From Eqs. ([T3|) to ([15|), Eq. ([23|), and Eq. ([24|) we obtain 

1 Q 1 

(25) 

It indicates that r2,2 = Ti is always in between Ti and T2 and smaller than the lesser of 4Ti/3 and 2T^. This 
prediction is similar to those obtained by others [23l. |24|. |65| and also agrees, within the experimental uncertainties, 
with the results of recent experiment. [14l | 

The relations between various characteristic times are summarized in TABLE |l| It is clearly shown that for a driven 
qubit the intrinsic decoherence time is independent of the driving field and is identical with the decoherence time of 
the qubit in free decay. In contrast, the relaxation and field-dependent decoherence times of the driven qubit strongly 
depend on the driving field and are equal to each other. 

TABLE L Relaxation and decoherence times of a qubit in free and driven decays. 
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I 1 


3 1 
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r2,2 


Ti 


2Ti 2T2 


4Ti 2T^ 



Decay type 


Relaxation 


Decoherence 




intrinsic field-dependent 


Free decay 
Driven decay 




T2,i = Ta T2,2 = Ti 



The relations between the relaxation and decoherence times of a qubit in free decay and driven oscillation can 
be used to extract the decoherence and dephasing times of the qubit by measuring its population evolution. In 
experiment, one would first measure the population evolution of the qubit in the free decay and in resonantly driven 
oscillation (Rabi oscillation) to obtain Ti and Ti. Then the decoherence times T2, T2.I7 and T2.2 can be evaluated 
using Eqs. ([^H)) to (I^S)) . Finally, the pure dephasing time of the qubit X^ can be calculated from Eq. ([T| 



V. MICROWAVE-DRIVEN SQUID FLUX QUBITS: ANALYTICAL VERSUS NUMERICAL RESULTS 

To demonstrate validity of the analytical expressions of relaxation and dccoherence times obtained in the preced- 
ing sections we apply them to calculate relaxation and decoherence times of a microwave-driven SQUID flux qubit 
inductively coupled to environment and compare the results with those obtained by numerically solving the master 
equation given by Eqs. H]) and ([2]). An rf SQUID consists of a superconducting loop of inductance L interrupted by a 
Josephson tunnel junction (J J) which, applying the resistively-shunted junction (RSJ) model ^66], is characterized by 
the critical current /c, shunt capacitance C, and shunt resistance R. The SQUID flux qubit is usually coupled to its 
control and readout circuits. A typical equivalent circuit is shown in FIG. 1 (a) along with its equivalent admittance 
Y{ijj) in FIG. 1 (b). In this simplified external circuit the left part is the SQUID and the right part supplies external 
flux to the SQUID qubit. For a superconducting device such as a SQUID, the thermal bath is the external circuit 
coupled to the device. Thus the external circuit is the dominant source of dissipation for the SQUID qubit [l^. In this 
section, we analyze the relaxation and decoherence of the SQUID flux qubit due to coupling to the external circuit. 

A. Hamiltonian of a microwave-driven SQUID flux qubit 

The Hamiltonian of a flux-biased rf SQUID with total magnetic flux <& enclosed in the loop can be written as [22.l67| 

H^{x)^^^V{x), (26) 

where, m — C$q is the mass of "flux" particle, $0 = /i/2e is the flux quantum, e is the elementary charge, x = $/<I>o 
is the canonical coordinate of "flux" particle, p = —ihd/dx is the canonical momentum conjugate to cc, and V{x) is 
the potential energy given by 

V (x) = -mujic (^ - ^ef - Ej cos {2ttx) . (27) 

Here, Ej = hlc/2e = ?7ia;|p/3^/47r^ is the Josephson coupling energy, /3^ — 2ttLIc/^o is the potential shape parame- 
ter, cjlc = 1/VLC is the characteristic frequency of the SQUID, and x^ = ^e/^o is the normalized external fluxes 
from the external circuit. 

To perform gate operations, a microwave pulse is applied to the SQUID qubit. If the interaction between the 
microwave and external circuit is negligible the Hamiltonian of the microwave-driven SQUID qubit coupled to the 
external circuit is given by 

H{x,t)^HQ{x)+HF{x,t)+Hi, (28) 

where, Hp (x, t) is the interaction Hamiltonian between the SQUID qubit and microwave and Hi is the interaction 
Hamiltonian between the SQUID qubit and the external circuit (thermal bath). 

If (f) (t) is the normalized flux coupled to the SQUID from the microwave then Hp {x, t) is given by [67| 

Hp{x,t)^'^^^^^[^ + 2{x~x,)]. (29) 

Hereafter 4> is taken to be 

m = cj)^ cos {u;^t) , (30) 

where, (p^ and ui^ are the field strength and frequency of the microwave, respectively. From Eqs. ([29|) and ([30|) one 
has 

muPjf^d),^ IX12I 
£= ^^J"^ ' , (31) 



and 



h 

For a SQUID qubit, 5 7^ if the tunneling distance jxn — X22I 7^ 



™^ic<^M (^11 -^22) 



B. Spectral density of the external circuit 

As has been addressed, in the case of weak dampmg, the effect of thermal bath on the superconducting system can 
be described by the spectral density. [62, l63| For the SQUID flux qubit considered here A = — l/$o and the spectral 
density J (w) is given by [6l|, ^M 



J (lu) ~ hLuYi^{uj) 



I + coth 



2kBT 



(33) 



where, kg is the Boltzmann constant, T is the temperature of the thermal bath, and Yfi{LL!) is the real part of the 
frequency-dependent admittance Y[uj) of the external circuit. [63| For the external circuit of FIG. 1 one has (see 
Appendix C for details) 

Yuu^ /f;:\ y (34) 

uj^ + Go ((jj) 

where, Fo (w) and Go (w) are given by Eqs. (|C8|) to (|C10|) in Appendix C. For the sake of concreteness, henceforth, 
unless otherwise specified, the SQUID flux qubit has the parameters oi L = 1.0 nH, G = 15 fF, Lg = 1.0 nH, Ge = 10 
pF, i?e = 10 17, i?o = 1.0 kil, [ij^ — 1.4, Xe = 0.499, and M = 3.0 pH. The temperature of the external circuit is taken 
to be T = 0.1 K. 

In FIG. 2 we plot the spectral density J [uj) versus frequency u for the SQUID qubit with the above parameters 
at T = 0.1 K. It is shown that the spectral density J (w) reaches the maximum at a; ~ 0.04iu;lc and approaches a 
constant value at low frequency, in particular J (0) 7^ 0. From Eq. (|18p T^p is finite and thus the external circuit will 
induce dissipation as well as dephasing. 

In FIG. 3, we show the spectral density J {uj) dX uj — lolc versus temperature (the solid curve) together with the 
quantum spectral density at T = (the dashed line) and the classical spectral density at high temperature limit 
T 3> hujLc/kB (the dashed dotted line). The quantum spectral density is independent of temperature while the 
classical spectral density is a linear function of temperature, which agree with those in the literature. [60| Note that 
for T > 1.5 K the spectral density is approximated well by the classical spectral density, while T < 0.3 K the spectral 
density approaches the quantum spectral density. 

For the SQUID flux qubit, Yr (w) is an even function of w, which is also the case for most of superconducting 
systems. Substituting Eq. ([55]) into Eqs. pB|) and (US]), we obtain 

T-i ^ '^hiu^i \x,2\'Yn{uj2i)coth (^^^ , (35) 



and 



T-^^^kBT{xii-X22fYHiO). (36) 



Thus the relaxation rate is dominated by the circuit's admittance at the transition frequency UJ21 while the dephasing 
rate by the admittance at w = 0. Furthermore, the dephasing rate is proportional to the temperature of thermal bath. 
Hence at the low temperature the dominating source of decoherence is relaxation while at the high temperature the 
main source of decoherence is dephasing. These results agree with those obtained by others. .30 , 56 , 62] 

C. Relaxation and decoherence times 

To numerically calculate the relaxation and decoherence times of the microwave-driven SQUID qubit, we need to 
compute evolution of population and coherence of the qubit. For this purpose we first calculate the eigenenergy i?„ 
and eigenstate \m) of the SQUID qubit by numerically solving the time- independent Schrodinger equation with Hq {x). 
[22ll67l| We then calculate the transition frequency comn, damping rate matrix element Rmn,m'n', and matrix element 
of interaction Hamiltonian i?,^„ between the SQUID qubit and microwave. Finally we calculate the population and 
coherence by numerically solving the master equation given by Eqs. ^ and ^ using the split-operator method [68| 
for the algebra equation with non-symmetric matrix. j61i] 



1. Free decay of the SQUID qubit 



In general, a SQUID qubit is a multilevel system. 67[ However, when it is driven by a weak resonant microwave 
field, the leakage to non-computational states is negligible [5J,l5a] and the SQUID qubit can be very well approximated 
by a TLS consisting of the lowest two levels |1) and |2) as the computational states. 

If the TLS is initially in an eigenstate (the ground or excited state) the coherence of the system will remain zero 
in free decay. To extract the relaxation and decoherence times of the SQUID qubit from numerically simulated time 
evolution in free decay we assume the initial state of the qubit is a superposition state with pn (0) = Pi2 (0) = 
P21 (0) = P22 (0) = 0-5- From Eqs. @ to (O, in this case uq = 1 and vq = wq = 0. Using this initial state one 
can calculate population and coherence of the SQUID qubit in free decay by numerically solving the master equation 
given by Eqs. IJ) and Q. 

In FIG. 4 (a) and (b) the solid lines are the numerical results of the evolution of population inversion {P22 ~ Pw) 
and squared modulus of coherence \pi2\ , respectively. Since the coherence is usually a complex and fast oscillating 
quantity, we use \pi2\ instead of Pi2 to estimate the decoherence time. It is shown that both the population inversion 
and the squared modulus of coherence undergo simple exponential decays. 

To evaluate the relaxation and decoherence times, we fit the numerically obtained {P22 ~Pii) and \pi2\ with 
exponential functions 

P22 -Pii =yi +2ie~'^^S (37) 



and 



Pi2?^y2 + Z2e-^"^\ (38) 



respectively. The results of least-square fitting are plotted in FIG. 4 (a) and (b) with dashed lines, where ti = 102.43 
ns and T2 = 188.27 ns (and other fitting parameters are j/i — —0.61198, zi — 0.61198, j/2 = 0, and Z2 = 0.25). For 
the weak damping considered here 6q w 0. According to Eqs. (fTO|) to (fT2|l. for the qubit in free decay, the relaxation 
time is Ti = ti and the decoherence time is T2 = T2- Using Eq. ([TS]), the dephasing time T^p = 2324.83 ns for this 
special case. Note that since in this particular case T^ ^ 2Ti the decoherence time is limited by the relaxation time. 
For comparison, we calculate the relaxation, decoherence, and dephasing times using the analytical expressions 
given by Eqs. (HH) to dUl), Eq. ^, and Eq. ^. The results are Ti = 102.43 ns, T2 = 188.27 ns, and T^ = 2324.83 
ns. These results are exactly the same as the numerical results, demonstrating the analytical and numerical methods 
give identical relaxation and decoherence times for the SQUID qubit in free decay within the relative error of the 
least-square fitting. 

2. Resonantly driven SQUID qubit - Rabi oscillation 

To illustrate the effect of driving field on relaxation and decoherence we consider a SQUID qubit being in the ground 
state initially driven by a resonant microwave field. For the ground state, pn (0) = 1 and P22 (0) — P21 (0) — P12 (0) = 
and thus uq = vq ~ and wq = 1 from Eqs. ([4]) to (l6|). In this case, from Eq. (|A7P the intrinsic decoherence terms of 
Eqs. (fTO|) and (P0|) vanish. With this initial state one can calculate the population and coherence of the microwave- 
driven SQUID flux qubit by numerically solving the master equation. 

In FIG. 5 (a) and (b) we plot with the solid lines the evolution of population difference {pn — P22) and squared 
modulus of coherence \pi2\ , respectively, for the SQUID flux qubit driven by a microwave with if) — 1.0 x 10""* 
and LOf^ = UJ21 — 7.22 x IO^^wlc, where lolc = 2.582 x 10^^ rad/s. The damping strength is 77 = 8.62 x IO^^ulc 
for the SQUID qubit considered here while the field strength corresponding to the critically-damped regime is </)„,, — 
4.94 X 10~^ from Eq. ([?T|) . Thus the SQUID qubit is in the underdamped regime. As a consequence, both the 
{Pii — P22) and \pi2\ undergo damped Rabi oscillations with Rabi frequency f2 = 1.74 x 10~^ ^LCi as shown in FIG. 
5 (a) and (b). 

To extract the relaxation and decoherence times of the driven qubit, we fit the numerically calculated {pn — P22) 
and \pi2\ to exponentially damped oscillating functions 

Pii - P22 =yi+zi sin (nt + ip^) e"*/^i , (39) 

and 

IPial^ = y2+Z2Sm{nt + ip2)e-*^^'^ 

+I3 sin^ {nt + (^2) e-2*/^22 ^ (40) 



respectively. The results of best fit are shown in FIG. 5 (a) and (b) with dashed lines, where ti = 132.68 ns 
and ?22 = 132.65 ns (and other fitting parameters are yi = 0.00033, zi = 0.99976, ip-^ = 1.56276, y2 = 0.00042, 
Z2 — 0.00665, Z3 = 0.24941, and tf2 — —0.01306). Since ^ <^ uj^ the system is in the weak field regime. According to 
Eqs. p9|) to ((2T|) . for the resonantly driven qubit, the relaxation time is Ti = ri and the field-dependent decoherence 

time is T2_2 = ^22. 

For comparison, we have calculated the relaxation and field-dependent decoherence times using the analytical 
expressions given by Eq. ([25|l and Eqs. (fT6|) to (flSjl . The results, Ti = r2_2 = 132.68 ns, are in excellent agreement 
with the numerical results. We have also calculated the relaxation and decoherence times of the underdamped qubit 
driven by microwaves of different field strength. The results are given in TABLE [Til It is shown that compared to free 
decay the effect of resonant microwave field is to make the relaxation time longer but the decoherence time shorter. It 
is also shown that the relaxation and field-dependent decoherence times obtained from the numerical calculation are 
essentially identical and independent of the field strength, which accord with the analytical results. Note that as the 
microwave field becomes stronger the numerically obtained decoherence and relaxation times begin to deviate from 
the analytical results due to strong field effects. 

TABLE II: Numerical results of relaxation and decoherence times (ns) of the SQUID qubit in free decay and driven decay. 



Field strength 


Relaxation time 


Decoherence time 





102.43 


188.27 


1 X IQ-*^ 


132.68 


132.68 


5 X 10"^ 


132.68 


132.68 


1 X 10"^ 


132.68 


132.67 


5 X 10"^ 


132.68 


132.67 


1 X 10"* 


132.68 


132.65 


2 X 10"'' 


132.70 


132.49 



VI. CONCLUSION 

In summary, by analytically solving the Bloch equation of a resonantly driven dissipative TLS under the RWA, the 
dynamical behavior of a driven qubit is systemically investigated. It is shown that the driving field has significant 
effect on the relaxation and decoherence of a qubit. For a resonantly driven qubit, the population and coherence 
undergo more complicated damped oscillations that in general have more than one exponential decay terms. Multiple 
relaxation and decoherence times are thus required to completely characterize the time evolution of the driven qubit. 
It is found that the decoherence of a driven qubit can always be decomposed into the intrinsic and field-dependent 
ones. The intrinsic and field-dependent decoherence times equal to the decoherence time of the qubit in free decay 
and relaxation time of the driven qubit, respectively. The relaxation time Ti and field-dependent decoherence time 
T2,2 of the driven qubit are always in between Ti and T2 and smaller than the lesser of 4Ti/3 and 2T,^. The analytical 
expressions for calculation of the relaxation and decoherence times of the driven qubit are derived. These analytical 
expressions have been used to study relaxation and decoherence of the microwave-driven SQUID qubit. The results 
are in excellent agreement with those obtained by numerically solving the master equation, confirming the validity of 
the analytical expressions. The relations between the relaxation and decoherence times of a qubit in free decay and 
driven damped oscillation are obtained. These relations can be used to extract the decoherence and dephasing times 
of a qubit by measuring its population evolution with free decay and resonantly driven Rabi oscillation experiments. 
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APPENDIX A: THE SOLUTION OF BLOCH EQUATION OF A DISSIPATIVE TLS IN FREE DECAY 

In the absence of ac driving fields (free decay), £==(5 = 0, and the field-dependent terms in the Bloch equitation ([7]) 
to ([9]) vanish. Usually, the coherence components u and v contain fast oscillating factors. Hence the first term on the 



right-hand side of Eq. ([9]) oscillates rapidly. In the rotating-reference which rotates with the coherence components, 
the last four terms on the right-hand side of Eq. ([7]) also oscillate fast. After dropping out all the fast oscillating 
terms under the RWA, the Bloch equation is approximated to 



du 

H 
dv 

H 
dw 

Itt 



{Rl2,12+ Rl2,2l)u + UJ2lV, (Al) 

-CJ21M+ (i?12,12 - -Rl2,2l)w, (A2) 

— (-^22,11 + -^11,22) W + (i?ll,22 — -^22, 11) ■ (A3) 



-Rll,22 


— R22SI 


-^22,11 


+ ^11,22 


y^ll- 


-^12,21' 


1 /^21 


+ ^12,21 



Obviously, there is no coupling between the population and coherence in this special case. Because of this the 
population and coherence evolve independently. 

If the initial value of the Bloch vector is (mq, vq, wq), the solution of Eqs. (jAip to (|A3p is then given by Eqs. (fTU]) 
to (|12p . The parameters and constants in the solution are given by 

(A4) 

(A5) 



2tan-^/ — '7"- -;, (A6) 

y CJ21 — -Ki2,2i ^ 

A^^, (A7) 

cos 00 

B=^, (A8) 

cos 00 

and 

D — Wo — Woo- (A9) 

APPENDIX B: SOLUTIONS OF BLOCH EQUATION OF A RESONANTLY DRIVEN DISSIPATIVE TLS 

For a driven dissipative TLS, e 7^ and also S ^ in general, as will be shown in section V. In this case, we assume 
that the trial solution of the Bloch equation ([7]) to (O still has the form of Eqs. PU)) to (fT^ but with the constants 
A, B, and D replaced by time-dependent variables /i, v, and A. Substituting the trial solution into Eqs. ([7]) to Q 
we obtain three equations with respect to the variables fi, v, and A. In the case of the resonant driving field with 
m^ — ccj, using the RWA, these equations are simplified to the form 

f - »■ <-' 

4^ = i/e cos (00)62"*, (23) 

dt 

where, 77 = (7 — k) /2 is the damping strength. Eqs. (|B2p and (jBSp indicate that due to the driving field the 
population and coherence are coupled to each other. 1^6^] Due to this coupling they no longer evolve independently. 

The form of the solution of Eqs. (|Bip to (jBSp depends on the relative value of e to 77. According to this relative 
value the dynamics of the driven dissipative TLS can be categorized into three regimes: underdamped regime when 
e > rj, critically-damped regime when s = rj, and overdamped regime when e < rj. For the underdamped regime with 
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e > r], the most important scheme for quantum computation, the solution of Bloch equation is given by Eqs. (fT9|) to 
()21l) . The parameters and constants in the solution are given by 

Bo = -^, (B4) 

cos 6^0 



\wo-Kx\ /-, 1 / vo + ex\ ,„,. 
cos(6'o) V ^ V Wo- KxJ 

Do = «:x, (B6) 

Di=Bicos9o, (B7) 

01=02+^3, (B8) 

^2= tan- ^i-^o-nx) ^^^^ 

evo - rjwo + {'qn + e^) x 

03 = tan^i -, (BIO) 



and 



7Wo 



e"^ + Kj' 



(Bll) 



For the overdamped regime with e < rj, the solution of Bloch equation can be obtained by replacing 51 with i\lr\^ — e^ 
in Eqs. P^ to (PT|) . In this case, the population difference and coherence components evolve with more than one 
exponential decaying terms and the coherence components are composed of intrinsic and field-dependent decohcrcnces. 
For the critically-damped regime with e = rj, the solution of Bloch equation can be obtained by setting fl = in Eqs. 
pop to (|2ip . In this case, the system undergoes a nonexponential decay owing to the nonexponential decay factor 
te~'"* in the population difference and coherence components. Note that the analytical solutions obtained here have 
been confirmed by the matrix exponent method. [69[ 

APPENDIX C: EQUIVALENT ADMITTANCE OF THE EXTERNAL CIRCUIT 

In this paper we assume that the circuit denoted by Ye{uj) in FIG. 1 (a) is an RC circuit with effective admittance 
Ye{uj) given by 

Y,{u;) = ^+ ^\. ^ . (CI) 

Ro Re + l/]UjCe 

The corresponding impedance Zg is given by 

Ye{uj) ]uCeff{uj) 

where, -Re// and Cg// are the effective resistance and capacity given by 

n / N Ro+LO C^ReRo{Re+ Ro) /^o^ 

1 +UJ^Ci [Re + Ro) 



and 



Ce//(^) — C'e 



1 I Re + Ro 



j^RlC^. \ R 
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(C4) 



respectively. They are functions of w. When Rq ^ oo, Reff = Re and Ceff = Ce 
For the external circuit shown in FIG. 1 (a), the circuit equations are 



Vi = jujLh^juMh, (C5) 

- -JLoAIh + (jUjLe + Reff + . ], ) h- (C6) 

The equivalent impedance Z (w) is therefore calculated from 

^M^^^,^^+ .,,/";,.^ . (C7) 

h J^Le + Reff + l/jLOCeff 

The equivalent admittance Y (uj) can be calculated from the equivalent impedance Z (lo) by Y (lu) = 1/Z (lu). The 
real part of the equivalent admittance, Y/j {u), is given by Eq. (|M)) . where, 

Fo{.)^^^, (C8) 

and 

c; = Af2 - LLe. (CIO) 

Eq. ([M|) shows that Yr (cj) is an even function of uj. It is also the case for most of superconducting qubits. 
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Figure Captions 

FIG. 1 (a) A simplified external circuit of the SQUID flux qubit. The left part is the SQUID qubit and the right 
part is the circuit used to supply external flux to the SQUID qubit. /i and Vi are the current and voltage of the 
SQUID qubit's circuit, I2 is the current of the external circuit, Le is the inductance of the superconducting loop of 
the external circuit, M is the mutual inductance of the SQUID qubit and external circuit, and Ye{uj) is the effective 
admittance of other devices in the external circuit, (b) The equivalent admittance Y (lu) of the external circuit. 

FIG. 2 (Onhne color) J{uj) vs. w for the SQUID flux qubit at T = 0.1 K. 

FIG. 3 (Online color) J(w) vs. T for the SQUID flux qubit at w = ujlc- The dashed and dashed dotted lines are 
the spectral densities at low temperature limit (T -^ 0) and high temperature limit (T -^ 00), respectively. 

FIG. 4 (Online color) Evolution of (a) the population inversion and (b) the absolutely squared coherence of the 
SQUID flux qubit in free decay. The solid and dashed lines are the numerical and fitting results, respectively. 

FIG. 5 (Online color) Evolution of (a) the population difference and (b) the absolutely squared coherence of the 
microwave-driven SQUID flux qubit. The solid and dashed lines are the numerical and fitting results, respectively. 
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